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ABSTRACT 

We compute the one loop vacuum polarization from massless, minimally 
coupled scalar QED in a locally de Sitter background. Gauge invariance 
is maintained through the use of dimensional regularization, whereas con- 
formal invariance is explicitly broken by the scalar kinetic term as well as 
through the conformal anomaly. A fully renormalized result is obtained. The 
one loop corrections to the linearized, effective field equations do not vanish 
when evaluated on-shell. In fact the on-shell one loop correction depends 
quadratically on the inflationary scale factor, similar to a photon mass. The 
contribution from the conformal anomaly is insignificant by comparison. 
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1 Introduction 



One of many exciting and potentially testable consequences of inflationary 
cosmology is a mechanism for generating the primordial cosmic magnetic 
fields which may have served as the seeds for the currently observed galactic 
field of about 10~^ Gauss. The idea [|l], |] is that the inflationary production 
of light, minimally coupled, charged scalars — such as the Higgs — resulted in 
the photon acquiring a plasma mass of about ~ eH, where H ~ 10^^ GeV 
is the inflationary Hubble parameter. Of course this would suppress the 
creation of photons during inflation, but it would vastly amplify the zero 
point energy of the super-horizon modes, 

—fiuj > —■ 

2 2 



-hu — > -^/Fe-2^* + m2 , (1) 



where k = 27r/A is the co-moving wave number. After the end of inflation 
the charged plasma dissipates — either by annihilation or through being 
redshifted into insignificance. If this happens quickly enough the enormous 
zero point energies are shed as coherent ensembles of very long wave length 
photons which would manifest as magnetic fields on super-horizon scales. 

To be more quantitative let us model the spacetime geometry during in- 
flation as locally de Sitter. We can express the invariant element conveniently 
either in co-moving or conformal coordinates, 

ds"^ = —dt^ + e^^^'dx ■ dx = a? —d-rf + dx ■ dx . (2) 

The conformal factor and the transformation which relate the two coordinate 
systems are. 

Gravity is a non-dynamical background. The dynamical variables are the 
vector potential v4^(x) and a complex scalar (f){x). Their Lagrangian is, 

C = -^F.^F.^g'^^g"'^^ - {d, - ieA,W{d, + ieA,)<Pg^^ y/^ ■ (4) 

One way of understanding the mass generation mechanism is by appealing 
to the result of Vilenkin and Ford for the coincidence limit of a free scalar in 
Bunch-Davies vacuum [131, 



^*(-)^(-)|^Le=(^) {UV + f^t}. (5) 
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Here "UV" stands for an ultraviolet divergent constant. If we infer an ap- 
proximate action for the photons by replacing the various 0*0 terms by the 
finite part of their vacuum expectation values, the result seems to be a time- 
dependent photon mass, 



Although very suggestive, the analysis of the preceding paragraph is not 
really consistent quantum field theory. The kinematic properties of a particle 
are encoded in its 2-point IPI (one-particle-irreducible) function. The 2- 
photon IPI function is known as "the vacuum polarization" and the one loop 
contributions to it are depicted in Figures 1-3. Making the replacement. 



corresponds to including only the diagram of Fig. 1. The other two graphs 
are the same order (e^) in perturbation theory and there seems to be no 
good reason for ignoring them. The diagram of Fig. 2 is required to make 
the vacuum polarization gauge invariant. And the graph of Fig. 3 is needed 
to absorb the ultraviolet divergence. 

The purpose of this paper is to compute all three diagrams in a consistent 
regularization and to demonstrate that they induce corrections to the photon 
wave function very like those of a photon mass. In Section 2 we review the 
familiar results from fiat space. This serves as a useful introduction to using 
dimensional regularization in position space and establishes a crucial corre- 
spondence limit for checking the accuracy of our subsequent work. Section 
3 gives the scalar propagator in Z)-dimensional de Sitter space. In Section 4 
we first reduce the vacuum polarization to manifestly transverse form, then 
we renormalize it. In Section 5 we take the result on-shell to demonstrate 
that the tree order wave functions receive one loop corrections like those of 
a photon mass. We discuss the result in Section 6, giving special emphasis 
to the important issues which are still open. A fuller discussion of what our 
result means physically can be found in another work [Q. 





e'A^A, 



2 Vacuum polarization in flat space 
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Figure 1: One loop contribution to the vacuum polarization from the 4-point 
(seagull) interaction. 




Figure 2: One loop contribution to the vacuum polarization from two 3-point 
interactions. 



8Z 

Figure 3: Photon field strength renormalization counterterm. 
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2.1 Momentum space treatment 

Although working in a curved background typically requires a position space 
treatment, it is desirable to begin our discussion of the vacuum polarization 
with the traditional momentum space treatment. The contribution from the 
4-point interaction is depicted in Fig. 1. It actually vanishes for massless 
scalar QED in dimensional regularization, 

(Fig.l) =-2.eV7p^jr— = »■ (8) 
Fig. 2 shows the contribution from two 3-point interactions, 

= (g^ "'"rgr'^' (p')*"(pV--PV). (10) 

Since we have the result for any spacetime dimension D we can note that 
its vanishing on-shell (i.e., at = in the transverse direction) for D = 4 
indicates that the photon stays massless in that dimension. The fact that the 
vacuum polarization does not vanish on-shell for D = 2 indicates that the 
massless scalars form a massive photon bound state in two dimensions. This 
is a scalar version of the Schwinger model ^ in which the photon develops 
a mass of ~ e. 

It is straightforward to extract the divergence as D = 4 — e approaches 4, 

(Fig. 2)"^ = -^{pY'-p''p')-^+ finite . (11) 
It is of course removed by the field strength renormalization of Fig. 3, 

(Fig. 3)^'' = -Z(5Z(r7^V - p^'p'') . (12) 
One sets AZ = — e^/x""/ {24:n'^e) plus any convenient finite term. 

2.2 Position space treatment 

Let us now repeat the exercise in position space. The propagator is, 

^A(x;x') = ^^-^^-^, (13) 

An 2 (Ax ) 2 
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where we define Ax^ to be, 



Ax^ =\\x-x'f - {\t~t'\-i5Y . (14) 

We will use the same figures for the position space diagrams. Of course the 
first diagram still gives zero, 



(Fig. l)'''' = -2ie^r]^"'iA{x; x)6^{x - x') = . 



(15) 



Note the general rule for using dimensional regularization in position space: 
the parameter e is assumed to lie in the range for which terms such as Ax"^"*"^ 
vanish at coincidence. For example, to derive (p!5| ) we assume e > 2, even 
though we shall eventually take e to zero. 

The position space versions of the other two diagrams are. 



(Fig. 2^ = ^e^ri^Pif" \dpiA{x; x')d'JA{x; x') - iA{x; x')dpd'JA{x; x'] 



(Fig. 3) 



i6Z{7]^''d' ■ d - d'^d'')6^ 



[X — X 



(16) 
(17) 



Note that these particular diagrams involve no integrations when written in 
position space. At two loop order there are integrations. 

The first step in evaluating (|^) is to substitute the propagator and take 
the derivatives, 



2e' 



r(i- 












(2-6 






'Ax'-^' 


rT2(l 


-f)(2 





i2-e) 



,Ax>'Ax'' 



Ax 



8~2e 



+ (2-e)(4-6) 



Ax 



6-2e 



+ 



'a^ 

2Ax'^Ax^' 



Ax* 



-2e 



(18) 
(19) 



We anticipate the notation of curved space in denoting the vacuum polariza- 
tion as a bi-tensor function, 
(/i) transforms according to t 



■'-'■1+2 



x; x'). This means that the first index 
le tangent space at the first argument (x) and 
the same relative relation exists between the second index (z/) transforms and 
the second argument (x'). 
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The next step is reaching manifestly transverse form. This is done by 
writing inverse powers of Ax as derivatives of lower powers. For example, 
one can easily derive the identity, 



1 



Ax 



4-2e 



-(4 - 2e) 



Ax 



6-2e 



+ (6 - 2e) 



Ax^Ax^ 



Ax^ 



-2e 



By combining with rj^^ times the trace it follows that. 



Ax^ 



-2e 



+ 



2Ax''Ax'^ 



Ax 



8-2e 



1 



2(2-e)(3-6) 



1 



Ax 



4-2e 



(20) 



(21) 



Substitution into (|l^) gives a manifestly transverse form. 



■■■^1+2 



x; X 



167r4 



3-e 



Ax' 



4-2e 



(22) 



The next step is extracting the ultraviolet divergence. This typically 
comes from a term of the form 1 / Ax^"^*^ through the identity. 



-2e(l 



Ax 



2-2e 



Ax 



4-2e 



(23) 



The reader will note that the various derivatives have so far failed to induce 
any delta functions. That was because there were either too few derivatives or 
because the power of Ax was wrong. Getting a delta function in dimensional 
regularization requires two derivatives acting on precisely the power Ax^~ , 



92 



Ax 



2-e 



(24) 



Combining (|23|) and (24) allows us to extract the ultraviolet divergence in a 
form which can be canceled by a local counterterm. 



Ax 



4-2e 



2e(l - e) Ax 



2-2e ' 



26(1-6) 

9^ rin(/i^Ax 



4 



Ax^-^^ Ax^-^ 

2\1 



6(i-e)r(i-f: 



(25) 

-5^(x-x'),(26) 



Ax^ 



e(l-e)r(l-f) 



5^(x - x') . 



(27) 
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As usual, dimensional regularization has resulted in a scale /x. 

Substituting (^TD reveals the divergence structure we already found in 
momentum space, 



1927r4 



87r2 V;,2 



e)(3-e) 



7]^"'d' ■ d - d'f'd'' 6^{x - x') , (28) 



247r2 



^""0' ■ d - d'^'d'' 



6^{x — x') — h finite . 
e 



(29) 



The cleanest subtraction in position space is to absorb the entire local term 
by choosing 5Z to be. 



5Z 



7T 



T(l- 



With this choice the fully renormalized vacuum polarization becomes, 

^\n{fi'^Ax^y 



(30) 



[x: X 



Ax' 



(31) 



2.3 Going on-shell in position space 

The position space derivation we have just completed was simpler than its 
momentum space cousin because no integrations had to be performed. They 
occur when one goes on-shell. What we are really checking, when we take the 
vacuum polarization on-shell, is whether or not there are quantum correc- 
tions to the linearized wave functions. One might define these as the matrix 
element of the operator Afj_[x) between the vacuum and a plane wave photon 
state. The equations obeyed by this matrix element come from varying and 
linearizing the effective action, 

T[A] = -^ld^xF,,F'^'' 

+ ^ J (fx J d^x'A^ix) [^W] (x; x')A,{x') + O(A^) . (32) 
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The associated field equations are 
6T[A] 



d^F^^'ix) + / (fx 



{x-x')A^{x') + 0{A^) = Q . (33) 



The general classical solution comes from superposing plane waves of the 
form, 

Al{x) = e,{k)e'- , (34) 

where = \\k\\ and the Lorentz gauge polarization vectors obey eo = = k-e. 
To see if there are quantum corrections one merely expands the solution in 
powers of h, 

A,{x) = A'^ix) + Alix) + . . . . (35) 

and then segregates all terms of the same order. Potential one loop correc- 
tions are determined by the equation, 



(fx^ 



{x-x')Al{x'). 



(36) 



We therefore conclude that the necessary and sufficient condition for there to 
he one loop corrections to the photon wave function is a nonvanishing integral 
for the one loop vacuum polarization against a classical plane wave solution. 



3 Scalar propagator in de Sitter D-space 

The behavior of free, massless and minimally coupled scalars on a locally de 
Sitter background has been investigated extensively 0, ||, 0; §]• Among the 
curious properties of these particles are the absence of normalizable, de Sitter 
invariant states and the appearance of acausal infrared singularities when 
the Bunch-Davies vacuum is used with infinite spatial surfaces [0, |[|. To 
regulate this infrared problem we work on the manifold T^^^ x R, with the 
spatial coordinates in the finite range, —H~^/2 < x^ < /2. Although the 
actual propagator is a mode sum on this manifold, the small possible varia- 
tion in conformal coordinates renders the first term of the Euler-Maclaurin 
formula — just the integral — an excellent approximation. So the finite 
spatial range of T^~^ serves merely to cut off what would have been a log- 
arithmic infrared divergence on R^~^. In = 3 + 1 spacetime dimensions 
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the result is 

iA{x; x ) 



D=4, 



H \ ^ f 1 1 1 



where the modified de Sitter length function has the definition,^] 



y{x; x') = a{ri)a{ri')H — — {\ri — r]' \ — iS) 



(37) 



(38) 



Neglecting the higher order Euler-Maclaurin terms does not prevent (^) 
from solving the correct differential equation. The higher terms also drop out 
of quite complicated, nonlinear relations such as the Ward identity for the 
one loop graviton self-energy . We shall therefore regard the technique as 
valid and confine ourselves to finding the appropriate generalization of (p7|) 
to D spacetime dimensions. 

We seek a function of y{x] x') and the two conformal factors which obeys, 



^^/^- Aa^-2(^) A,A(x; x') = z5^(x - x') . 



(39) 



When the kinetic operator acts on a function of just y{x] x') one finds, 

= H\'>{^){{Ay - y^)f"(y) + D{2 - y)f(y) - iiS(n " '/)/'(!/)^}. (40) 

The only symmetric function of 0(77) and 0(77') which can give the same 
prefactor of a^{ri) is a constant times the same logarithm that appears in 
). The D-dimensional propagator must therefore take the form. 



iA(x;x') = f[y{x]x')^ + 61n(^a(?7)a(ry')j , 



(41) 



^What is termed "the de Sitter length function" in the literature is, 

x') = 1 — y{x] x') . 

The geodesic length from to x''^, £{x; x'), is related to y{x; x') as follows, 



y{x; x') = sin^ (^^H£{x; x')^ 
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where the function f{y) obeys, 

H'a^iv){iAy~y')ny)+D{2 
-Ai5{r^-il')f{y)5-h{D- 



y)f\y) 



(42) 



The deha function (for 5^0) obviously descends from a factor of ?/^~~ 
Series solution of the equation then generates an infinite sum of higher pow- 
ers. Defining D = 4 — e and normalizing correctly gives, 



47r^ 



-rr 1 




1 



E 



r(3 + n- f) 



4n+l 



(43) 

This series solves (|4^ ) for 6 = 0, but it does not reduce to (|37|) for e = 0. The 
n = term of the sum is not even finite in this limit! The resolution to both 
problems is a series of strictly nonnegative integer powers of which cancels 
the divergence and the unwanted terms. This series obeys the homogeneous 
equation up to a constant which is canceled by the b{D — 1) term. 



iA(x; x') 




1 



r(3 + r2- 



+ 



n + f (n + l)!r(2-f 




The great advantage of this regularization is that it preserves general co- 
ordinate invariance (once 5 is taken to zero). One might think that the prop- 
agator is unwieldy but this is not so in practice. For example, this formalism 
has recently been used to compute and renormalize all two loop contributions 
to the stress-energy tensor of a real scalar with a 0^ self-interaction |TlI]. The 



really cumbersome part of (|4^) is the infinite sum on the second line. But 
these terms all vanish at coincidence {y{x\ x) = 0) and they vanish for all 
y{x; x') at D = 4. So one need only retain them when they multiply some- 
thing else that diverges like 1/e. Note also that one need never worry about 
large y{x; x') on account of causality. 

All valid regularizations must reproduce the result of Vilenkin and Ford 
that the coincidence limit of the propagator contains a finite term which 
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grows like ln(a) = Ht To check this note that y{x;x) = at coinci- 
dence. When a variable vanishes like this in dimensional regularization one 
must always assume e to be large enough that the variable is raised to only 
nonnegative powers. We therefore find, 

hm,A(.;.')= (|:)'(^)"*{^r(3- j) + lr(3-.)ln(a(„))} . 

(45) 

Note that is exact for arbitrary e and indeed reduces to in the limit 
when e approaches zero. 



4 Vacuum polarization in de Sitter 



The diagrams which contribute to the vacuum polarization in de Sitter back- 
ground are drawn exactly as those of fiat space, so we shall use the same 
figures. Of course there are now some factors of the de Sitter metric! We 
express this in conformal coordinates and adopt the usual convention that 
indices are raised and lowered by the Lorentz metric. 



(Fig. l) 



flU 



-2ie^^-g{x)g^"'{x)tA{x; x)6^{x - x') , (46) 
-2ie^a°-^ri'"'iA{x; x)S^{x - x') . (47) 



(Fig. 2)'" = 2e^^-g{x)g^P{x)^-g{x')g''''{x') 

dpiA{x; x')d'^iA{x; x) — iA{x; x')dp&^iA{x; x] 
'd^'tAix; x')d"'iA{x; x') - tA{x; x')d^'d"'iA{x; x') 



X 



(Fig. 3) 



-t5Zd,(^ [g^'^gp'^ - g^'^g-'P] d'Xix - x' 
-i5Z [7]>"'d' ■ d - d'^'d"] a^~^5°{x - x') . 



(48) 
• (49) 

(50) 
(51) 



In these and subsequent expressions we define a = a{r]) and a' = a{r]'). Note 
as well that the scalar propagator is the de Sitter one (PI). 

It is convenient to subsume the complicated, e-dependent constants which 
appear in the propagator (^), 

(52) 



iA{x;x') = a^'y(y{x; x')^ + l3\n(^a{ri)a{ri')^^ . 
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Comparison with (0) reveals, 



a = 



27lJ 



r(2- 



2^ ' 



1 r(3 - 6) 
4r(2-f^ 



(53) 



and, 



7(y) ^ , . . + (1 - i) 



1 

X) 

+ E 



2?/ ^ 

°° lT{3 + n-e) y" 



n=l 



1 r(3 + n 



n r(2 + n - e) 4"+^ n + f (n + l)!r(2 - f ) 4"+i 



In this notation derivatives of the propagator ( ^2]) can be written as. 



dy 



d,d'^tA{x;x') = a{Y{y) 



„, N dy dy 



+ i{y)- 



d'^y 



dxP dx'"' ' ' dxPdx'" J 
We can therefore express the portion of (^8]) within the brackets as. 



(54) 

(55) 
(56) 
(57) 



a' 



dpiA{x; x')d'^iA{x; x') — iA{x; x')dp&^iA{x] x' 

dy dy 



7"-7"(7 + /51n(aa') 



dxP dx'"' 



7 



7 + /3 In (aa' 



d^y 



a6' 



dy , , dy 



dxPdx"" 
+ (3^H^aa'6\6\ . (58) 



''dx"' dxP " 

It is straightforward to differentiate the de Sitter length function ([38| ) 

dy 



dxP 
dy 
dx'" 
d^y 
dxPdx'" 



a{r])H [y6^^ + 2a{r]')HAxp + 2iHa{r]')sgn{r] - r]')6^p6\ , (59) 
a{ri')H \y5°^ - 2a{ri)HAx^ - 2iHa{ri)sgn{ri - rOb'^A , (60) 



aa'H^ 



y6\6\ - 2a{Ti)8\HAx, + 2a{ri')HAxp6\ - 2r]p 



(7 -^v'//" p--— ^cr I -^v'/ /^^ "-^P" CT -'/per 

2m(r/)a(r/')/^^|r/ - VI V-*^ - ^^^iv - v')S\S°j] ■ (61) 
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Gauge invariance requires taking 6 to zero. Only the final order 6 term can 
contribute in this limit, and then only when it multiplies 7', 



— Aiaa H'^6{ri — rj') Mm^y' (jj{x; x')^6 



5^ (a; - x') 



aa 



D-2 



(62) 



The net result is to subtract off the purely temporal components of (^71) . We 
accordingly combine this term with (Fig. Ij to form, 



+ ^|34yMa)}5^(^-^')- (63) 



A bar over a tensor indicates its zero components have been removed. 

The left-over portion of (^Fig. 2^^^ will be known as i '^Ilg (x;x'). 

Expression (|58|) seems complicated due to the infinite sum in the definition 
of 7(2/). However, we need only retain terms which survive as e 0, 



1" - 7" 



7 + /?ln(aa') 



1 



1 - f y'-' 



+ 



,3-e 



[2- 



e \2 



'2-f\ r(3 - e) In(aa') , \Yi{H^/\x^^ 



7 



2e y 
7 + /3 In(aa') 



3-1 



y3 2 



+ 4?' 



1 1 ^(l-e)(2-|) 1 



1 - f y3-e 



26(1-1) 1/2- 



1 1 r(3 - e) In(aa') ^ ln(i72Aa;2) 



V 2e ; 7/2-2 4r(2- f) 7/2-2 ■ 4?/ 
After taking e to zero in other non-divergent terms the result is 



(66) 



[x- x') = 2e^a^a^-'a'^-' I Aa^a'^H'^Ax^Ax" 



1 - f y'-' 



(2-1)2 1 



r(3 - e) In(aa') ln(i72Ax2^ 



2e y 
-2aa'H^r]^"' 



3-4 



V 4 ; r(2 - f ) y3-f 
1 1 (l-e)(2-|) 



1 - f 



2e(l - 



,2-e 



+ 



4^2 



1 

4y2 



3-e 
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1 r(3 - e) In(aa') \n{H^A 



+ 4r(2-f) 
\n{H^Ax^) + 1 



X 



Ay 



+aa'H'^5%5 



\n{H^Ax^ 



(67) 



2^,2 t7-2e-4 



7]- 



1 /2- 



e VI - 







Ax^ 






i) 


aa' 










2 









6-2e 



Ax 



(1 - e)r]^"' - (4 - 2e) 



Ax^Ax' 



Ax' 



Ax' 



,4-2e 



1 r(3-e) 

2r(2-f; 



In(aa') 



X 



(4-e: 



Ax'^Ax'^ 
Ax^ 



Ax 



4-e 



2 /2 

-a a 



-a-^HAx^d^o + a'-^6%HAx'' - H^Ax^Ax' 



\n{H^Ax^) + 1 
A? 



16^ 



+ 



2 /2 

a a 



ln(/f2Ax^ 
A?~ 



■ (68) 



The next step is reducing to manifestly transverse form. The first term in 
8]) is exactly the same as the flat space result ([19D and its reduction gives 
HI), as before. The second term is reduced with the identity, 



(1 - e)ri^"' - (4 - 2e) 



Ax'' Ax' 
Ax^ 



Ax^-2^ " 



2-2e" ' ~ ' ~ ' Ax^-^' ■ 

„e-4 



(69) 



The fact that the third term goes like Ax^ means that its reduction pro- 
duces a local term. 



(4 



Ax^Ax' 



1 



Ax' 



4-e 



2 -e 



2m^-f 
f{2 



Ax^-' 
14 



r7^^(5^(x - x') . (70) 



This local term completely cancels 
1 and 2 is therefore, 



1+2 



[x; x') 
1 

rjT]' 

+ 



TrT^fl - 



2(3 - e) 



). The sum of all terms from Figures 
1 



Ax 



4-2e 



\n{H^Ax^) + l 



Ax' 



+ r]- 



■fiu 



\n{H^Ax^ 
'2rfrfAx' 



T]Ax^5\ - id^Ax" - Ax^Ax' 



\n{H^Ax^) + 1 
rfri''^Ax 



(71) 



Note that only the first term — the one that survives in the flat space limit 
— still requires regularization. 

Although expression (^) is transverse this is not yet quite manifest owing 
to the factor oil/ still standing to the left of the derivatives in the second 
term. After some tedious tensor algebra one finds, 



1+2 



e2 r7rT2(l-f) 



2(3 



Ax^-^' 



■\yi{H^Ax') + 1 
r]r]'Ax'^ 



+ 



|ln2(if2Aa;2) + l\n{H^Ax' 



(72) 



It remains only to extract the divergence and subtract it with Fig. 3. Al- 
though the divergence is exactly the same as in flat space (p8|) the countert- 
erm is not, owing to the factor of a^~^ in (51). The incomplete cancellation 
gives rise to a finite factor of ln(a) times a delta function. Of course this is 
the conformal anomaly |jT2|, |13[ . The renormalized vacuum polarization is. 



{x; x') 



in-" 



/ln(/i2A 



X 



V 2AAx' 



+ 



\n{H^Ax^ 



TjTj'Ax 



— ln(a)(5^(x — x) 
3 



+ 



TjM-V' ■ V - d'^'d" 



■\n\H^Ax^) + l\n{H^Ax^ 



.(73) 



Note that it is completely integrable and gauge invariant, and that we have 
everywhere taken D = 4. 
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5 Going on-shell 



Taking the vacuum polarization "on-shell" in a locally de Sitter background 
is complicated by the fact that the "in" vacuum is not equal to the "out" 
vacuum. This is obvious from the fact that there is particle creation. The 
photon wave function is therefore not the matrix element of An{x) between 
a 1-photon in-state and the out-vacuum. It is rather the matrix element of 
A^{x) between a 1-photon state and the Bunch-Davies vacuum, both pre- 
pared at t = 0. The field equations obeyed by this matrix element are given 
by varying the Schwinger-Keldysh effective action [l^, [1^ . 

The rules for computing in the Schwinger-Keldysh formalism are simple. 
The diagrams have the same topology as those of Feynman but the endpoints 
of lines bear either a "+" or a "— " polarity. All external lines are +, whereas 
vertices can be either all + or all — . The + vertices are the same as those 
of the standard Feynman rules; the — vertices are conjugated. There are 

-|— 1-, H — , — h and propagators. All of them are the same function 

( ^11) of the appropriate version of the modified de Sitter length function 
y±±{x;x') = a{ri)a{ri')H^Ax^^, where we define, 

Ax1^{x;x') = \\x — xW"^ — {\t] — rj'l — id)"^ , (74) 

Axl^{x;x') = \\x - x'f - {t] -r]' + i6f , (75) 

Ax'^_^(x;x') = \\x — x'W'^ — {r] — 7]' — i5Y , (76) 

Axl^{x;x') = \\x-x'f - {\T]-T]'\+i6f . (77) 

^A_^_^{x■, x') and iA_^._{x■, x') are equal for rj' > r], hence the ++ and H — 
contributions cancel whenever rj' > rj. When the x'^ and x''^ are spacelike 
related, the real part of y{x; x') is positive; when they are timelike, the real 
part of y{x;x') is negative. Therefore meromorphic functions of y^^{x;x') 
agree for spacelike separation for 6—^0. That is why the ++ and H — 
contributions cancel when x'^' strays outside the past lightcone of x'^. Inside 
the past lightcone the ++ and H — propagators are conjugate. 

The Schwinger-Keldysh effective action involves background fields A^(x) 
for the -|- lines and A~{x) for the — lines, 

r[A+; A-] = -J I F+ - F-,F-4^7^^^7-v^ 




+ d^xd^x'\A^{x)[^'Ul^ {x;x')At{x') + A^{x) ^'Ul_ {x; x')A- {x') 
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[x] 



{x-,x')A-{x')\ + 0{A').{7S) 



The various ± permutations of the vacuum polarization are all the same 
function ( |73D with the appropriate ± permutation ([7^[77D substituted for 
"Ax^" The Schwinger-Keldysh field equations are obtained by varying the 
effective action with respect to either background and then equating the + 
and — fields, 



5V[A+;A- 



per 



+ 



-^11. ^/J 



rfV { ['^n:;^] (x; x') + ['^n:;.] {x- x')]a^{x') + 0{a^) 



[x \X] 







[X\ X 



(79) 



Note that we have exploited the relation 

Because electromagnetism is conforma 
field equations are the same (in conformal coordinates) as those of flat space. 
The general classical solution is therefore a superposition of plane waves, 



ly invariant for D = A the order 



ik-x 



(80) 



where k!^ = \\k\\ and the Lorentz gauge polarization vectors obey = = k-e. 
As in flat space we check for quantum corrections by expanding the solution 
in powers of h, 

A,ix) = Alix) + Alix) + . . . . (81) 

and then segregating all terms of the same order in the field equations. Po- 
tential one loop corrections are determined by the equation, 



\2 iiv 



Alix) = - j rfV { [^Yil^^ {x- x') + ['^n';_] (x; x')]Al{x') . 



(82) 

^This may seem surprising because there are no H — or — h seagull graphs or coun- 

{x\ x') or its conjugate is 



terterms, so the only diagram topology contributing to 
Fig. 2. The apparent paradox is resolved by noting that the mixed permutations also 
fail to produce local terms coming from partial integration. Owing to the absence of the 
temporal absolute value in Aa;^_ we have, 



Ax 



,4-2e 



9^ 

"2e(l-e) 



Axi 



Axi 



' 4 



\nin^Axl_) 
Axi 
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The temporal integration in this equation begins at the initial time of r( = 
—H~^. Its upper limit is irrelevant as long as it comes later that rj because 
the cancellation between ++ and H — contributions eliminates contributions 
from any point x'^ which is outside the past lightcone of x^. 

From (^) we see that the vacuum polarization can be expressed as the 
sum of four distinct terms, 



\^iC J tXj ^ 



le 



F{Ax 



2 , , G{Ax±±] 



+ 



=F— ln(a)5±±(5 {x - x') 



O O 



rjr]' 
H{Ax. 



2 1 



(83) 



Note that the conformal anomaly only contributes for the ++ and cases. 

The functions of Ax\j^ in the other three terms are, 



F{Ax^) 
G{Ax^) 



d' 



\n{iJ,^Ax'^ 



2AAx'^ 
^\n{H^Ax^) + 1 
A? 



192 ^ 96 . 

— \n\H'Ax^) + -\n(H^Ax^] 
16 ^ ^ 8 ^ 



K{Ax^) = -\n''{H^Ax^) + -\n{H^Ax^) 
8 2 



(84) 

(85) 
(86) 
and 



It is useful to make the additional definitions, F(Ax^) = d'^f{Ax^ 
G{Ax^)=d^g{Ax^). 

Of course it is easy to compute the conformal anomaly's contribution to 
the right hand side of (|82D . Since this term is local we can actually derive it 
for an arbitrary vector potential and then specialize to the classical solution, 



C^(x) = 



247r2 

^2 



T]'"'d' ■d-d'^'d" \n{a)5\x-x')\A,{x'), (87) 



247r2 
24^ 



5^(ln(a)F'^^(x; 



Hka{ri)e^'{k)e 



ik-x 



(88) 
(89) 



C^lx) completely dominates the classical term by virtue of the factor of «(?/), 
but it is much less significant than the a'^{r]) associated with a true photon 
mass. Interestingly, the conformal contribution is also purely dispersive. 
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Since the anomaly contribution (^) is local, it is instructive to combine 
it with the classical action 5*0 into an effective action, 



X 



1 + 



247r- 



ln(^)U^ 



X — X 



^A,{x'). (90) 



More generally, for a nonabelian gauge theory G coupled to Nf Dirac fermions 
and Ng complex scalars in representation r, we have []T^ [T^ |n 



5'o + 5S'a 



X 



1 + + ^^/C*!^) - ^^C2{G)) In ( 



ao 



V^^v'^F^uFp,, (91) 



where g is the gauge coupling constant, ag = g'^/4:Tc, C2{G) is the quadratic 
Casimir of the adjoint representation of G, which can be defined in terms 
of the generators of the adjoint representation T*^ and the group structure 
constants /'^"'^ as tr[T'^T''] = -f^^^^f^^ = -C2{G)S''\ and C(r) is defined in 
terms of the generators t"(r) of representation r as tr[t"(r)t^(r)] = C{r)6"''^. 
For example, for SU(N), G2 = N, and for fermions and scalars in the fun- 
damental representation of SU{N), G{r = fundamental) = 1/2. The effect 
of the anomaly (|9l| ) on the photon dynamics is quite moderate Since 
the photon mass contribution is parametrically much larger, we expect the 
photon dynamics to be affected much more by the photon mass |]l|. 

Our strategy for evaluating the other three contributions is to exploit the 
fact that the range of x''^ is independent of x'^ to first pull the derivatives 
outside the integration. We then combine the ++ and H — terms to obtain an 
integrand which is only nonzero for x''^ within the past lightcone of x'^. The 
final step is a sometimes lengthy series of asymptotic expansions under the 
assumption that the mode under study went super-horizon long ago (hence 
k -C Ha{ri) = after a long period of inflation (hence H <^ k). We 

shall organize these expansions in terms of two dimensionless parameters. 



y 



-krj 



and 



w = 



k 
H 



(92) 



The physically interesting region is 
crossing one would have y ~ 1 and w 



<y<til<^w<^a. At horizon 
^ a. Evolution of the physical scales 
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in de Sitter inflation is illustrated in figure 4. The physical wave length 
Aphys = ^na/k = 27!- /Hy grows with time, while the Hubble radius remains 
constant. 




n 



'0 'X 

Figure 4: Evolution of the physical scales in de Sitter inflation. Horizon crossing 
occurs at rjx- 



Although it is the least important of the four terms, an exact result can 
be obtained for the "F" contribution to the right hand side of (p^). We begin 
by reflecting the x'^ derivatives {d'^F{Ax'^) = —d^F^Ax"^)) and then pulling 
all derivatives outside the integration. 



F^ix) 



It 



ri^^'d' ■ d - d'^d"] \f{AxI^) - F{Ax 



>eue 



ik-x' 



ik-x' 



(93) 
(94) 



The difference between the ++ and H — contributions vanishes for x'^ outside 
the past lightcone of x^ , 



I'K f r 

/(A4+)-/(Ax+_) = -{ln[/.^(Ar^ 



-\\e{Ai])e{Ari-\x\) . (95) 



We next change variables to Ax^ = — x'^ and perform the angular inte- 
grations, 

-)2 _ . 



F^ix) 



le 
8^ 



dArje 



ikAr) 
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x47r / drr — — - — 
io kr 48 



{ln[At2(A?72-r2)] -1} . (96) 



The radial integration involves a combination of special functions ^ (/c A77) 
that we will meet again, 

j drr sin(A;r)|ln /i^(A7y^ — r^) — l| 

= A77^y" rfa;a;sin(A;A77x){21n(/xA77) - 1 + ln(l - x^)} , (97) 

sin(A;A?7) - A;A?7COs(A;A77)] [21n(/xA?7) - l] + Ari'^^{kAri) .(98) 



1 

fc2 



Our definition for ^ is," 



e(«) 



= / (ixx sin(Q;x) ln(l — x ) , 

Jo 

2 1 
= — sin(Q;) r-[cos(Q;) + a sin(Q;)] 



si(2a) + I 



H — -[sin(Q;) — q;cos(q;)] 



ci(2a)-7-ln(|) 



(99) 



(100) 



Here 7 ~ .577 is Euler's constant. The following small and large a expansions 
are sometimes useful, 



5-5ln(2) 



a + 0{a^) , 



COs(q!) + - 

a. OL 



7COs(q;) — — sin(Q!) 



+ 



ln(Q;) 



(101) 
(102) 



Combining the previous results gives F^{x) = e^d^ — e ■ dd^ d^e^'^I{rj) 
where. 



dArje 



ikArj 



^We use the following notation for the sine and cosine integrals: 

sin(i) 



si(x) = 
ci(x) = - 

J X 



t 



TT 
"2 



dt- 







t 



°° , cos(t) , , , 

dt ^-^ = 7 + ln(a;) 



dt 



cos{t) — 1 
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X { [sin(A;A77) - kAr]cos{kAr])^ [21n(/xAr;) - l] + k'^Ar]'^^{kAr])} . (103) 

The only x dependence in this expression resides in the outer exponential fac- 
tor. We can act with the spatial derivatives and then commute the temporal 
derivatives through the exponential to obtain, 



2\3ik-x ; 



(52 - 2tkdor'i{v) . 



(104) 

(105) 



After commutation the only 77 dependence to the right of the derivatives 
resides in the upper limit of the At] integration. So one of the 9o's acts to 
undo the integral and we are left with the messy task of taking the remaining 
derivatives. The final result is, 



kH 



247r2 1 - a-i 



k 



1 — a 



_gi|-(l-o-i) I ^IJ-^ik-x 



(106) 



Because it lacks positive powers of the scale factor, F^{x) is far weaker 
than C'^(a;). In fact it is only nonzero because the process begins at t = 
(77 = —H~^), rather than at t = —00 [r] = —00). This is as it should be since 
the "F" contributions to the vacuum polarization are the same as those of flat 
space, which vanish when exactly on-shell. We get a nonzero result because 
beginning at a finite time precludes one from going precisely on-shell. This 
is an important check on the accuracy and consistency of the process. 

The "K" contributions are reduced in almost the same way. They have 
zero temporal component by definition. Their spatial components are. 



K'(x) 



Stt^ 



(jij V' • V - &'d^ 



H{Axl^) - H{Axl_) 



>eje*'^"',(107) 



Stt^ 



e'V' - e • Vd"] I d^x'^[H{Axl^) - H{Axl_)] , (108) 



i-ri+H-^ pikA-q 

JO [At] — riY 



x47r 



2Sill(^"^) ^TT 



drr 



kr 



{ln[ii'2(A7y2-r2)] +2} , 



(109) 
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iV72 



ik-x 



ri+H- 



dAr] 



{Arj — rf)' 



sin(A;A?7) - kAr] cos{kAr]) 2 \n{HAr]) + 2 + k^Ar]^^{kAr]) \ , (110) 



p2 tj2 pw—y p 



-a e e 



Att"^ Jo {z + uY 

X {[sin(z) -2cos(z)][21n(z) -21nH + 2] + z^i{z)] . (Ill) 

(Recall that we define y = —krj -C 1 and w = k/H ^ 1.) The next step 
would be making an asymptotic expansion of the integral but this would be 
wasted effort because almost the same integral occurs with the opposite sign 
in the "G" contributions. 



The contributions from G ( |85D to the right hand side of ( p2|) are the 
most difficult to evaluate owing to the combination of temporal derivatives 
and the factor of 1/rjrj'. Because of this it is desirable to act with some of 
the derivatives and exploit the transversality of the polarization vector at an 
earlier stage than with the other terms, 



le 



1 



G{Axl^) - G{Axt 



X / d xe 



4^1 „—ikAxQl 



lg2^ik-x 
V 

'g{Axl 



Tjr]' 



giAxl_) 



ik-x' 



(112) 



rj' 



(113) 



X / d xe 



+ Vi-S'odo + ik")] {do - 2tk)do 
'giAxl^) - g{Axl_) 



i^i —ik-Ax a' 



(114) 



At this point it is best to partially integrate the d'^. Owing to the cancellation 
between g{Ax\j^) and g{Ax\_) outside the past light cone, only the lower, 
temporal surface term can survive. 



4^1 „—ik-Ax Qf 



d X e 



g{Axt 
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-ik-Ax 



7]' 



(115) 



r]'=-H- 



Since k'^k^y = we see that only the i/ = component survives. 
The spatial integral can be evaluated by familiar techniques, 



d-'x'e 



g{Axl^) - g{Axl_ 

sin(A;Ar7) — kArj coslkArj) 



(116) 



^ikArj 



X 



= -p-S(/cA?7,w) . 



2 \n{kAT]) - 2 ln(u;) + 1 + k^Ar]^^{kAr]) } , (117) 



(118) 



The function S(a;,w) has the following asymptotic expansions for small and 
large x respectively, 



K , J 3 



In(^) -1^ +0(a;^ln(x)) , 



w 6 



(119) 



xe ■ 



(21n(w) -ln(2x) - (1 - 7)) cos(x) sin(x)j + 0(\n{xw)^ .{120) 



Substituting the spatial integral into ( |115| ), and inserting the result into 
(ITT4D gives. 



G^'ix) 



^2^Meifc-^(l _ ^5^) (2 - idy)dy I ^ dz 



87r2 



87r2 



z + y 
w 



(121) 



2^f^^ik-x^2-idy-2ydy + iyd^y) / dz 

<J 



{z + yf 



"8^ 



'e^e''-^{l - tdy - 2ydy + lyd^^ (2 - id, 
24 



, r^^-^'^^ (122) 



(Recall that we define y = —krj <^ 1 and w = k/H ^ 1.) Now note that the 
first term on the top line almost exactly cancels i^'^(x), 



87r2 



'^-y , 'Eiz.w) 
az- 



{z + yf 



az 



(z + yy 



smu — 2;cosu 



dz— 

z^ 



smu 



z cosu 



e^'^-"{ln(w) + 0(1)} . 



(123) 

+ 0(y)| , (124) 
(125) 



This turns out to be the magnitude of the leading order contribution for 
super-horizon modes late during inflation. 

Another leading order contribution comes from the second term on the 
top line of (|122| ). To save space we suppress the prefactor of ^^a^e'^e*'^'^, 



id,, 



az- 



{z + yf 



:^(w -y,w) + 2i 



{z + yf 



(126) 



The flrst term on the left hand side is of order ln{w)/w but the integral can 
make a leading order contribution. Since the integrand converges for y = 
and w oo, the desired leading term comes from the explicit factor of \n{w) 
in the function S(z,w), 



2i 



{Z + |/)3 

— 4iln(w) / dz— 



sm[z) — zcos[z 

oo 



—4i\n(w) X — - 

^ ^ 4^2 

21n(ti;) + 0(1) . 



e^'' -l-2iz 



+ 0(1) 
+0(1) , 



(127) 

(128) 
(129) 



The flnal leading order contribution comes from the flrst term on the 
bottom line of ( |122| ). It is useful to flrst note the exact identity, 



{2 + id;,)E{x,w) 
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-x{ ci(2x) + ln(2x) + 1 - 7 - 2 \n{w) + i si(2x) + 2 } + e*^{sin(2;) 

si(2x) + |]}.(130) 



X [ci(2x) + ln(2x) + 8-7-2 ln(w) 
Setting X = w — y we conclude, 



t cos X 



(2 - tdy f^"^ J'"^' = InH + 0(1) . 



(131) 



Combining all of the leading terms we obtain the following result for the right 
hand side of Eq. (|82|) , 
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C^'ix) + Ff'ix) +G^'{x) + K^'ix) =a\r])——{\n{w) + 0{l)}A'^^{x) . (132) 

Note that the contributions from C'^(a;) and F^{x) are down by factors of 
w/a = fcphys/i^ -C 1 and {w/a)^, respectively. Equation ( [132| ) is consistent 
with a photon mass of, 



(133) 



^ 27r2 ""Vff. 

It is important to keep in mind that this result is perturbative. A full non- 



perturbative analysis of (|8^) is thus required in order to calculate reliably 
the photon mass. 



6 Discussion 

We have presented a long calculation in perturbative quantum field theory 
and it is worth commenting on the matter of reliability. Our result (^) for 
the vacuum polarization passes many important consistency checks. The first 
of these is gauge invariance. This is not trivial even in flat space. In a locally 
de Sitter background it requires a horrifying series of seemingly unrelated 
terms to combine into transverse projection operators. Yet they do. 

Another important accuracy check is that there are no new ultraviolet 
divergences. There should not be if the theory is to stay renormalizable 
because there are no new counterterms in de Sitter background. Related to 
this is the fact that our result has the correct fiat space limit. A final check 



is that the conformal anomaly term agrees with standard results |12, 13 . 



26 



The structure of the self energy ( |73| ) is very similar to that of the photon 
self energy in thermal QED. The self energy (in momentum space) is usually 
parametrized by two functions, Il^'^{k) = Pj^'^Ilt{k) + Pj^'^Ili{k) , where PJ^'^ = 
V^Y^iSij - kikj/P) and P^"" = t]^'' - k''k''/k'^ - P^" are the (spatially) 
transverse and 'longitudinal' (in fact time-like transverse) projectors. Our 
result ( [73D can be viewed as the space-time generalization of the transverse 
and 'logitudinal' vacuum polarizations, 

Ilt{x,x') = (9'-(9nW(x,x')+V'- Vn(2)(a;,x') 

Ui{x,x') = 9'-(9nW(x,x'), (134) 

where 11*^^) and 11*^^) are the transverse and spatially transverse contributions 
to ([7^) (c/. also Eq. ( |135D ). In thermal QED vacuum polarization modifies 
the photon dynamics. In the static limit 'longitudinal' photons are Debye- 
screened by the fermionic plasma, with the Debye mass, tjid = eT/^/3, 
while transverse photons are screened only dynamically. At high momenta 
k ^ T, the 'longitudinal' modes become unphysical, while the transverse 
photons propagate as massive particles with the thermally induced mass 
rriT = eT/v^. Based on the above mentioned similarities, we expect that a 
more detailed study of the vacuum polarization (|73D should reveal analogous 
physical effects on the photon dynamics in inflation. 

We are very confident about our result ( |73D for '^11'^ {x,x'). It has to 
be admitted that the process of going on-shell is much less well checked. 
This is also the most complicated part of the calculation. One important 
point is that (|106|) , the flat space contribution F^(x), vanishes in the limit 
that the initial time is taken to negative infinity. Although F^[x) makes the 
weakest of the various contributions, its reduction is quite similar to that of 
the crucial G^{x) and K^{x) contributions. So the fact that F^{x) obeys an 
important correspondence limit partially checks them as well. 

Our result is consistent with a photon mass of = ^\n{k/H). In- 
terestingly, this is precisely what follows from the Hartree-Fock estimate (|^) 
if one replaces the time dependent factor of ln(a) by its value at horizon 
crossing, \\i{k/H). It is premature to make too much of this coincidence. 
Although our 1-loop vacuum polarization is exact, all the work of taking 
it on-shell really demonstrates is that one loop corrections to the classical 
photon wave function become non-perturbatively large. 

To actually solve for the photon wave function and show that it ap- 
proaches that of a massive photon requires two extensions of the current 
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work. First, we must establish control over higher loop corrections. If the 
one loop correction becomes large then why does the two loop correction not 
give an even bigger effect? There is a curious parallel between what we must 
do and the problem that Schwinger faced for D = 2, massless QED in flat 
space @. Just as it was possible to use the one loop result in that context, 
so we believe it will prove to be in this case. 

To understand how this can be, note that we know quite a lot about the 
general structure of the vacuum polarization. As a consequence of gauge 
invariance, spatial translation invariance and spatial rotational invariance, it 
must have the form, 



{x;x') = [r]^"'d' ■d-d'^'d'']A{7],7]';Ax 



+ 



TjM^V'. V-a'V fi(r/,r/';Ax2) . (135) 



The functions A(j],ri'; Ax^^ and 3(^7], rj'; Ax^^ must be symmetric under in- 
terchange of 7] and 1]', and they must have the dimensions of inverse length 
to the fourth power. The full flat space result must reside in A(r],r)'; Ax'^^, 
but it will be negligible in de Sitter background. Terms that matter for de 
Sitter are those with factors of 1/t] and l/t]'. We believe quite strongly that 
A(ri, 1]'] Ax^^ can contain at most l/{riri') and that B(ri, rj'] Ax^^ can contain 
at most l/{rfri''^). 

Obviously each loop will contribute a factor of e^. The only really diffi- 
cult thing to guess is the number of logarithms. We believe — though less 
strongly — that the general result at i loops is A{ri, rj'; Ax^^ contains up to 

I logarithms whereas b(j],7]'; Ax^^ contains up to £ + 1. We further believe 
that in each case these logarithms can translate, after going on-shell, into up 
to one factor of ln{k/H) for each loop order. Since there is also a factor of 
a = e^/47r for each extra loop, this suggests that retaining only the one loop 
part would be reliable for modes which obey aln{k/H) < 1. That is, the 
one loop term would dominate the classical one by a factor of a^(?7), but the 
two loop correction to it would be down by a factor of aln{k/H). 

This seems a reasonable and probably provable conjecture. If it is true 
then we can essentially get the full photon wave function by solving the 
integral-differential equation obtained from just the classical term and the one 
loop vacuum polarization in (^). Which brings us to the second necessary 
extension of the current work: solving such an equation. Note that spatial 
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translation invariance implies we can still use spatial plane waves. But then 
the spatial integrations are identical to the ones we already did in Section 5! 
We can therefore reduce the problem to one of solving for the multiplicative 
function of 77. This is probably tractable analytically in some reasonable 
approximation. If not, it can certainly be done numerically. 

Another issue is the extent to which our process repeats the Higgs mech- 
anism of flat space. We do not have Poincare invariance during inflation 
but one's expectation is still that a massive photon has three polarizations. 
Since no new degrees of freedom have been granted to the vector potential 
one might suspect that the third polarization must come from the derivative 
of the scalar phase. It would be interesting to check this. 

Finally, we comment on the possibility of important stochastic effects. 
Although we have defined the photon wave function as the matrix element of 
Af^{x) between Bunch-Davies vacuum and a simultaneously prepared plane 
wave photon state, one should bear in mind that the quantum averaging 
implicit in matrix elements may give misleading results. The actual charge 
density induced by inflationary particle production is not smooth but rather 
stochastic [0, corresponding to a highly squeezed state. About one massless, 
charged scalar exists per Hubble volume, moving at the speed of light in a 
random direction. We do not expect that this distinction amounts to a 
significant difference for super-horizon photons, which are affected by many 
different Hubble volumes. This can and should be checked. 
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